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Abstract 

Low-temperature quantum oscillations of the dHvA amplitude are shown to 
persist far below the upper critical field of a strongly type-II superconductor, 
due to the gapless nature of the BCS quasiparticle spectrum in high fields. 
The dHvA amplitude in the superconducting state is smaller than its normal 
state counterpart by the factor ~ [max(T, r)/A]^, where T is the damping. 
This factor reflects the presence of a small gapless portion of the Fermi surface, 
surrounded by regions where the BCS gap is large. The agreement with recent 
experimental data on V^Si is very good. 
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There has been much interest lately in properties of high Tc superconductors (HTS) and 
other strongly type-II systems at high magnetic fields Recent reports clearly demonstrate 
the de Haas- van Alphen effect in the mixed state of A15 superconductors {VsSi and NbsSn), 
as well as in the layered superconductor 2H — NbSe2 0-3- In all cases it was found 
that quantum oscillations in the dHvA amplitude persist to a surprisingly low fraction (~ 
60%) of the upper critical field Hc2- This is surprising since, in the standard Abrikosov- 
Gor'kov theory, one expects an exponential suppression of dHvA amplitude due to the large 
superconducting gap A at the Fermi surface. Furthermore, it was found that the dHvA 
amplitude for the fixed value of magnetic field, H, behaves as a function of temperature, 
T, in the same way as in the normal state, except for the overall reduction in magnitude 
when the sample becomes superconducting. The presence of these dHvA oscillations in the 
mixed phase could, in principle, be due to a small portion of the sample remaining normal. 
However, heat capacity measurements have ruled out this possibility p|,||. 

In this Letter we show that these experimental results are a direct manifestation of a 
qualitatively new nature of the BCS quasiparticle spectrum at high fields. As shown in Refs. 
I^J^, at fields near Hc2, this spectrum is gapless at the set of points in the magnetic Brillouin 
zone (MBZ). These nodes in the gap refiect the center- of-mass motion of Cooper pairs in 
high magnetic field, in contrast to the familiar nodes of a zero-field unconventional [i.e., p- 
or d-wave] anisotropic superconductor which are due to the relative orbital motion. The 
gapless behavior persists to a relatively low fraction of Hc2 0, as long as A(T, H) remains 
smaller than or comparable to huc, where A(T, H) is the average BCS gap and Uc = eH/mc 
is the cyclotron frequency. The results of Corcoran et al. then follow from the presence of 
a small portion of the Fermi surface containing a coherent gapless band of quasiparticles, 
while the rest of it is gapped by a large A High magnetic fields, low temperatures and 
clean samples provide ideal conditions for the validity of this picture, with the coherent 
quasiparticle propagation extending over many unit cells of the vortex lattice. For example, 
the V^Si sample used in experiment satisfies well the last condition, with its electronic 
mean-free path, /q, being much longer than the separation between the vortices (given by the 
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magnetic length / = yc/2eiJ): Iq = 1450 A at T = 1.3K, while the intervortex separation 
varies from ~ 60A to ~ 45 A for fields between 10 and 18.5 Tesla. This situation should 
be contrasted with the one at low fields, where coherent propagation is suppressed and the 
low-lying quasiparticle excitations are localized in the cores of isolated vortices. 

At high magnetic fields, the electrons participating in superconducting pairing occupy a 
well-defined set of Landau levels (LLs). At low temperatures and high fields, where 3> 
T, r^^ (Fq = h/2T is the scattering rate) and huc ^ A, numerous 'quantum oscillation' effects 
arise as consequence of this Landau quantization In this regime the superconducting 
order parameter A(r) is well-described by the Abrikosov solution confined to the lowest 
LL for Cooper pairs [|10|- The BCS Hamiltonian is diagonalized in the representation of 



Ref. 1 11 1, where electronic states are labeled by the quasi-momentum q [q ± H] restricted 
to the first MBZ, momentum kz parallel to the field and the LL index n. Cooper pairs are 
formed from the electrons with the opposite momenta (q, k^) and spins belonging to the same 
(diagonal pairing) or different LLs (off-diagonal pairing) [§,0. Near Hc2, where A -C huc, 
the off-diagonal pairing can be treated as a perturbation and the BCS quasiparticle spectrum 
is obtained analytically |l5|J^JT0(]: 



Efc.,q,n = ±V^n(A;.)2 + |A„„(q)|2 (l) 

h'^k'^ 

Bnikz) = + fii^c{,n + 1/2) - /i, 

where /x is the chemical potential. The spectrum consist of branches (n^ = int[/i//ico'c] 
is the number of occupied Landau levels) in which the superconducting order parameter 
A„„(q) goes to zero at points qj, forming a 'vortex lattice' in the q-space. There is a strong 
linear dispersion around these points (nodes). A„„(q) for different n in (Q) behave similarly 
in the vicinity of and differ considerably only far away from these nodes. At lower fields, 
where A becomes comparable to TiuJc-, the off-diagonal pairing must be included on equal 
footing. The excitation spectrum is found numerically in Ref. [^] and still exhibits gapless 
behavior in some of the branches, while the gap starts opening up in the rest. Eventually, 
at yet lower fields, the LL structure is destroyed by large A. p The presence of these nodes 



and other regions with a very small gap at the Fermi surface, as A(T, H) increases from 
zero to ^ huJc, is the key aspect of our theory. 

We now consider the dHvA effect in such gapless high-field superconductor in the low 
temperature limit. The first step is to determine the oscillatory part of the thermodynamic 
potential, Qosc{H,T). Initially, we consider a homogeneous 3-D system [like V^Si, Nb^Sn 
and other A15 type-II superconductors] and then show how the results change for the layered 
systems like NbSe2. 

The quasiparticle contribution to the thermodynamic potential per unit volume is given 

by: 

^exc = E TTT- E E In [1 + exp(-/3E„(A;„ q))] + ^ [E^{h, q) - (2) 

y spin ^xJ-'yJ-'z „=o fc^q I ^ J 

where /3 = 1/T and En{kz,q) is the quasiparticle excitation energy (|l]). Using the standard 
Poisson resummation formula: 

POO roo 

V (f){n) = <f){0)/2 + / dn(j){n) + 2Re V / dn(j){n) exp{i27Tkn) (3) 
n=o ^0 f^^^Jo 

we can perform the sum over Landau level index n in (0) and get the oscillatory part of ^lexc 

as: 

1 ^ oo 27[k IT 



/ d^q I In 

Jmbz 



l + exp(-/3V^2 + |A(q)|2) 
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where q = ql is the quasi-momentum rescaled by magnetic length. 

In deriving (4) we have assumed a spherical Fermi surface. In real systems, with 
more complicated Fermi surfaces, cos {27c kfi/huc — 7r/4) is replaced by cos {hA/uj*m* — 7r/4), 
where A is the extremal cross-sectional area of the Fermi surface in the plane _L H and 
uj* = eH/m*c, m* being the effective mass associated with the orbit around A [0). We 
have also neglected the fact that A, the overall amplitude of A„„(q), oscillates as a function 
of H, and have assumed that the oscillatory piece of A is much smaller then its 'smooth' 
part. This is a good approximation if the number of occupied LLs, ric, is large [y/n^ ^ 1]. 



Furthermore, we have neglected different features in A„„(q) as a function of a LL index n 
since they are pronounced only away from the gapless points [nodes]. 

While the systems in question are rather clean 0], some impurities and defects are still 
present and give rise to a small exponential decay of dHvA amplitude even in the normal 
state. Therefore, we should include disorder in our calculations. This is accomplished by 
introducing the density of states per gapless branch of the spectrum (|l]) in presence of 
disorder as: 

W„(o)4*„g-5^|(^iLp . (5) 

where u = uj — S(u)) and Nn{0) is the normal density of states at the LL. In the clean limit 
(no disorder) N{uj)/Nn{0) ~ (c<j/A)^ for small energies, reflecting the presence of nodes 
(gapless points) in the quasiparticle spectrum. This strong energy dependence suggests that 
the effective scattering rate in the superconducting state will be modified relative to the 
normal scattering rate t~^. We take the self-energy S(cj) to be S((I;) = —iT, where F = h/2T, 
which results in Ndis{(^) / Nn{0) ~ (F/A)^ + (cij/A)^. This is a good approximation in the 
unitary limit (otherwise the scattering rate depends strongly on energy-we have investigated 
the behavior of the density of states in presence of various forms of static disorder and will 
present these results elsewhere [|l3l). In the end, F has to be determined self-consistently as 

It is known that in conventional superconductors, with a finite gap everywhere in the 
momentum space, the quasiparticle contribution to Qosc{H,T) is exponentially small at low 
temperatures because of the large value of A/T (unless the field is very close to Hc2, within 
few percent). On the other hand, in the gapless high- field superconductor described above, 
there are quasiparticle excitations with momenta q such that A(q) < T which give a large 
contribution to the thermodynamic potential (4). Therefore, we divide the MBZ in two 
different regions. The 'gapless' region (which we label Q), has A(q) = Aq for all q < qc, 
with (jc = Cmax(T/ A, F/A) being the radius of Q and C a constant of order unity. Outside 
of Q, where A(q) > max{T, F) for all q, we will aproximate A(q) with the average value of 
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the gap function A =< A(q) >av IHI- This is the 'gapped' region. In setting the boundaries 
of Q we have allowed the possibility of T < F, which is actually the case in Ref. |0. There, for 
example, T < 0.2 meV while the normal scattering rate for V^Si superconductor is Tq = 0.36 
meV in the normal state (we will show that the self-consistent F actually increases slightly 
from the normal state value). 

Inserting the density of states (H) in (4) and taking the derivative over H, we obtain the 
oscillatory part of the magnetization due to the qasiparticle excitations from the region Q 



as: 



eh 1 jj, 



T F 



-1 2 CO f -| Nfc 



^'-^^'^) = -;;^^(^ ^-^^(a'a) E^sin(^-./4) (6) 



•A' A 



1)'^ . ,27Tkn 



k-- 



2Ti^kT f T V \ 

X T/ sinh (— ) exp {-lixkV jhu^ + O [max( — , —)f 

nujc \ A A / 

The amplitude of dHvA oscillations A{H, T) of magnetization (6) for the first Fourier har- 
monic k=l is: 



T F 



2 



T 



AHH, T) oc 2 [Crna.(-, -)J H-^'^--^^^^-,^^ exp (-2.r/S..) (7) 

which is the same contribution as in the normal state reduced by the factor G = 
2[Cmax(T/A, F/A)]^ (for the moment we assume F ~ Fq). Factor 2 in (|^) comes from 
the presence of two gapless points in the First MBZ of the spectrum (p. The result (0) tells 
us that, in passing from the normal to the superconducting state, there should be a drop in 
the magnitude of dHvA oscillations, since now there is only a small fraction G ^ 1 of the 
Fermi surface where there are low-energy excitations. The size of Q is determined by both 
the total number of nodes at the Fermi surface and the areas in different branches where the 
BCS gap is very small although not necessarily zero. In presence of finite disorder both will 



broaden into gapless regions ||T^. At low temperatures, such that T < F, this broadening 
is quite small: only ~3% of Fermi surface in V^Si is still gapless at T = 1.3K by the time 
magnetic field drops from H = 18.5(~ Hc2) Tesla to = 10 Tesla. 

The quasiparticle excitations with momenta outside Q, where A > max{T, F) , contribute 
to the amplitude of dHvA oscillations as: 



A'-^{H, T) oc (1 - G)H-^/^^f{^^) exp - 2'KT/huo,) (8) 



where f{x) = 1 for x = and f{x) = y7rx/2 for x ^ 1 and we have assumed k = 1. This 
amphtude is much smaller then A^: it is already ~ 20 times smaller then (7) at if = 16.5 
Tesla in VsSi system and becomes negligible at lower fields. 

Figure 1 shows the field dependence of ln[A{H, T) sinh (27r'^T/hiUc)H^^'^T~^] [Dingle plot] 
for V^Si sample used in the experiment by Corcoran et al. [0, both in normal and super- 
conducting states. The normal state data (straight line) give the normal scattering rate 
of To = 0.36 meV. The dashed line in Figure 1 shows the theoretical plot obtained with 
A{H, T) = A^{H, T) + A^^^{H, T), where we have taken the scattering rate F equal to its 



normal value Fq. We have assumed A(i7) = A(0)y 1 — H/Hc2, which is a good approxi- 
mation for the range of fields used in the experiment. The comparison of this curve and 
the experimental points suggests that there is a small increase in the scattering rate in the 
mixed phase. Indeed, the self-consistent calculation of F gives ||15| : 



nH) = f-^ for F/A<1 (9) 

Using this value of F in (7) and (8) we obtain the full line in Figure 1. Our theoretical curve 
is in excellent agreement with the experimental data. We should note here that we use only 
a single fitting parameter, C ~ 1.012, which yields a reasonable value of ~ 3% for the size 
of Q. All other quantities are taken from Ref. 0. 

Finally, we address the problem of dHvA oscillations in layered superconductors, like 
2H — NbSe2 and HTS Y Ba2Cu^0T-s. The layered systems have been extensively studied 
ever since Graebner and Robbins ||^ reported first observation of dHvA oscillations in 2H — 
NbSe2 at magnetic fields below Hc2- In a layered system there is hopping, t, between 
the layers so that the electron dispersion along the field is given by t cos {kzd), where d is 
the interlayer separation. In this case the oscillatory part of the magnetization due to the 
quasiparticle excitations around gapless points is given by: 

MLiH^T) = -^y^T^-^TTT^ [cmax(|, I)] ' E ^ sin - ./4) (10) 



mc l^dihuj^flH^I^ 



k 



X T/sinh (-^^) exp {-2'KkT /hu:^) + O [[max{-, -)fj , 

where we have assumed that t > fiujc- Comparing (6) and (10) we see that magnetization in a 
layered system has the same form as in the homogeneous 3-D system up to the dimensionless 
factor {■\/miT^d)~^. The dHvA oscillations persist deep in the mixed state because of the 
gapless region Q on the Fermi surface. However, in the layered superconductor such as 
2H — NbSe2 there is a possibility of the charge density wave instability that can modify the 
shape of the Fermi surface in the certain parts of the MBZ, so that the results of the simple 
model presented above might not be quantitatively appropriate. 
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FIGURES 

FIG. 1. ln[^sinh(X)ijV2r-i] in V^Si as function of l/H, where X = 2iT'^T/nuc and 
A = A^{H,T) + A^~^ {H,T). Full circles represent the experimental data of Corcoran et al. 
[2]. The dashed line is a theoretical result obtained from Eqs. (7,8), with F set to Fq = 0.36 
meV, while the full line uses F evaluated self-consistently (9). The vertical dashed line indicates 
the normal-superconductor transition at i?c2 = 18.5 Tesla. We have used Aq = 2.6 meV and the 
effective mass of 1.6 irie and 1.7 irie in the normal and vortex state respectively, (from Ref. [3].) 
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